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AN  EXPLORATORY  STUDY  OF  A  THREE-DIMENSIONAL 
SHOCK  WAVE  BOUNDARY  LAYER  INTERACTION  AT  MACH  3 


by 

B.  Oskam,  I.  E.  Vas  and  S.  M.  Bogdonoff 
Princeton  University 
Princeton,  New  Jersey  08540,  U.S.A. 


SUMMARY 

An  exploratory  experimental  investigation  has  been  carried  out  on  the  three  dimensional  flow  fields 
caused  by  the  interaction  of  oblique  shock  waves  and  a  planar  turbulent  boundary  layer.  The  study  was 
performed  at  a  free  stream  Mach  number  of  2.95,  a  Reynolds  number  per  inch  of  1.6, x  10^  and  near  adiabatic 
wall  conditions.  The  interaction  was  studied  on  two  experimental  configurations  having  different  initial 
boundary  layer  thicknesses  (6  =  0.13  and  0.55  inches).  Both  surface  measurements  as  well  as  complete  flow 
field  surveys  were  performed. 

The  main  contributions  of  the  present  investigation  are  two  experimentally  derived  flow  field  models 
for  shock  generator  angles  of  4°  and  10°.  Based  upon  both  static  pressure  and  surface  flow  patterns,  as 
well  as  heat  transfer  data,  the  interaction  region  can  be  characterized  as  quasi -two-dimensional  along  the 
shock  direction  in  the  region  studied. 

A  critical  examination  of  the  occurance  of  ^ordinary"  flow  separation  and  its  character,  as  applied 
to  the  present  problem,  was  carried  out.  It  was  concluded  that  McCabe's  criterion,  as  used  by  the  pre¬ 
vious  investigators,  is  not  a  sufficient  condition  to  determine  the  onset  of  flow  separation. 
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1.  INTRODUCTION 


Since  the  early  studies  of  supersonic  flow,  there  has  been  continued  research  on  the  problems  of  the 
interaction  of  shock  waves  and  turbulent  boundary  layers.  These  studies  have  been  concerned  with  geo¬ 
metries  where  the  pressure  gradients  were  imposed  in  the  streamwise  direction,  resulting  in  planar,  two- 
dimensional  flows.  In  most  cases  of  practical  importance,  there  are  probably  significant  three-dimensional 
effects.  In  nominally  two-dimensional  flows  this  three-dimensionality  is  associated  with  end  or  comer 
restraints  and  involves  gradients  in  flow  properties  in  all  three  directions.  However,  there  is  a  large 
class  of  three-dimensional  flows  which  contain  one  particular  direction  along  which  flow  properties  vary 
only  slowly.  In  this  class  of  shock  wave/boundary  layer  interactions,  one  might  expect  a  quasi-two- 
dimensional  character  and  therefore  a  simplification  over  the  full  three-dimensional  problem. 


The  study  presented  herein  deals  with  such  a  quasi-two-dimensional  interaction  between  an  oblique 
shock  wave  and  the  turbulent  boundary  layer  on  a  surface  parallel  to  the  stream  but  perpendicular  to  the 
plane  of  the  shock.  Interactions  of  this  kind  are  found,  for  example,  on  the  sidewalls  of  rectangular 
supersonic  air  intakes  or  oblique  shock  diffusers.  The  geometry  produces  both  streamwise  as  well  as 
transverse  pressure  gradients  resulting  in  significant  cross  flows  in  the  boundary  layer.  These  complex 
three-dimensional  flows  have  been  the  subject  of  only  a  few  studies,  primarily  because  of  the  experimental 
difficulty  involved  in  the  full  resolution  of  the  flow  phenomena.  McCabe*  and  Lowrie^  have  studied  the 
problem  both  experimentally  and  theoretically.  McCabe  presented  an  approximate  inviscid  theory,  based  on 
vorticity  arguments,  which  relates  the  limiting  surface  streamline  angle  to  the  shock  conditions.  This 
theory  incorporates  some  fairly  sweeping  assumptions,  which  result  in  an  oversimplification  of  the  physi¬ 
cal  situation  and  therefore  severely  limits  the  understanding  that  can  be  gained  from  it.  A  more  compli¬ 
cated  and  also  more  realistic  treatment  has  been  given  by  Lowrie.  His  analysis  is  primarily  based  on  the 
argument  that  the  pressure  gradient  term  in  the  momentum  equation  is  much  larger  than  the  shear  stress 
term  and,  consequently,  the  problem  is  treated  as  inviscid  rotational  fluid  flow.  The  development  of 
cross  flows  within  the  boundary  layer  can  be  predicted  in  this  way,  but  only  in  an  initial  region  of  the 
interaction  where  the  pressure  gradient  is  indeed  the  dominant  effect.  It  should  be  noted  that  Lowrie 
assumed  constant  static  pressure  in  the  direction  normal  to  the  wall  in  both  his  theory  and  his  experi¬ 
ments.  He  reported  that  boundary  layer  profiles  that  showed  a  static  pressure  variation  of  more  than  5% 
were  discarded. 

Token3  and  Neuaann4  have  concentrated  on  surface  heat  transfer  as  it  is  most  important  for  direct 
practical  applications.  Although  McCabe,  Lowrie  and  Token  provide  a  large  amount  of  information  about  the 
three-dimensional  interaction  of  an  oblique  shock  wave  with  a  turbulent  boundary  layer  each  of  these 
studies  leaves  a  significant  part  of  the  phenomenon  unresolved. 


The  present  study  involves  the  same  general  problem  as  examined  by  the  previous  investigators  but 
with  a  somewhat  different  objective;  to  examine,  in  a  well-defined  experimental  set-up,  the  complete  flow 
field  of  the  interaction  and  to  directly  connect  this  flow  field  with  the  surface  phenomena.  The  goal  is 
to  provide  a  better  understanding  of  the  phenomena  by  considering  many  different  aspects  of  the  complete 
interaction  and  to  form  a  framework  for  future  experimental  and  theoretical  studies. 

Variable  strength  shock  waves,  at  a  Mach  number  of  2.95,  interacting  with  a  high  Reynolds  number 
turbulent  boundary  layer  were  the  general  framework  of  the  experimental  program.  The  shock  interaction 
studied  was  at  some  distance  from  the  shock  generator  to  try  to  determine  the  characteristics  of  the  inter¬ 
action  in  a  region  not  directly  influenced  by  the  streamwise  comer  formed  by  the  shock  generator  and  the 
test  surface.  Major  concentration  was  on  the  study  of  the  interaction  using  one  relatively  thick  turbu¬ 
lent  boundary  layer  (Model  1),  but  the  scope  of  the  study  was  extended  to  a  much  thinner  turbulent  boun¬ 
dary  layer  (Model  2)  to  examine  scaling  effects. 


Special  attention  was  paid  to  the  phenomenon  of  flow  separation  in  three  dimensions.  It  should  be 
noted  that  a  unique  characteristic  of  three-dimensional  flow  is  that  it  can  separate  from  the  surface 
without  the  mean  surface  shear  stress  becoming  tero.  This  type  of  flow  separation,  known  as  ordinary 
separation,  has  to  be  distinguished  from  singular  separation  whose  principal  characteristic  is  that  the 
mean  wall  shear  stress  is  tero  at  the  separation  point  (e.g.,  the  separation  in  planar,  two-dimensional 
flows  is  singular).  Previous  investigators  of  the  present  problem  have  indicated  the  occurance  of  ordi¬ 
nary  separation  along  a  line  defined  on  the  basis  of  surface  oil  flow  patterns.  The  present  investigation 
examines  this  ordinary  separation  on  the  basis  of  both  surface  phenomena  as  well  as  flow  field  measurements. 

2.  EXPERIMENTAL  PROGRAM 


2.1  Geometrical  Arrangements 

A  schematic  of  the  experimental  configuration,  showing  the  region  of  detailed  measurement,  is  pre¬ 
sented  in  Fig.  1.  The  oblique  shock  was  generated  by  a  variable  angle  shock  generator  which  spanned  the 
tunnel  between  top  and  bottom  wall.  As  noted,  the  region  of  study  extends  upstream  of  the  interaction  to 
establish  the  incoming  boundary  layer  characteristics.  The  test  region  does  not  extend  to  the  shock  gen¬ 
erator,  and  is  limited  in  extent  by  the  expansion  from  the  trailing  edge  of  the  shock  generator  and  the 
reflection  of  the  shock  from  the  wall  opposite  to  the  shock  generator.  In  any  such  experimental  set-up, 
extraneous  interference  effects  from  the  boundaries  of  the  region  of  study  can  be  expected  and  great  effort 
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was  made  to  evaluate  these  effects.  The  length  of  the  shock  generator  was  made  long  enough  so  that  the 
expansion  from  the  trailing  edge  did  not  interfere  with  the  region  of  study  and  a  detailed  examination  was 
made  of  the  region  of  influence  of  the  reflection  of  the  shock  wave  from  the  opposite  wall. 

One  of  the  key  parameters  of  the  study  is  the  effect  of  shock  strength  which  was  varied  by  changing 
the  deflection  angle  of  the  shock  generator.  Two  configurations  of  this  general  format  d  ig.  2)  were 
studied  in  detail;  their  main  difference  was  the  thickness  of  the  boundary  layer  coming  into  the  inter¬ 
action  region.  Most  of  the  studies  were  carried  out  on  Model  1,  because  the  relatively  thick  turbulent 
boundary  layer  allowed  high  resolution  measurements  to  be  made  throughout  the  interaction  region.  The 
configuration  of  Model  2  incorporates  a  much  smaller  incoming  boundary  layer  thickness  and  thus  covers  a 
larger  region  of  the  interaction  in  terms  of  boundary  layer  thickness  than  Model  1.  A  more  detailed  des¬ 
cription  of  the  experimental  geometries  is  given  in  Ref.  5. 

2.2  W ind  Tunne 1 


The  experimental  study  was  carried  out  in  the  Princeton  University  blowdown  wind  tunnel.  The  test 
section  has  a  cross  area  of  8  inches  by  8  inches.  The  test  conditions  for  the  present  tests  were  a  Mach 
number  of  2.95,  a  stagnation  pressure  of  100  psia,  an  average  stagnation  temperature  of  472°R  and  a  tunnel 
wall  temperature  of  10  to  20%  above  the  adiabatic  wall  temperature.  These  test  conditions  resulted  in  a 
Reynolds  number  per  inch  of  1.6  x  10^  and  a  maximum  test  time  of  6  minutes.  Details  of  the  basic  facility 
are  given  in  Ref.  6. 

Previous  measurements  had  indicated  that  the  boundary  layer  flow  over  the  bottom  and  top  walls  were 
equivalent  and  planar,  two-dimensional.  The  bottom  wall  boundary  layer  was  used  as  the  test  layer  for  the 
Model  1  study.  Windows  positioned  in  the  top  and  bottom  walls  were  used  either  as  observation  or  instru¬ 
mentation  ports.  In  general,  pressure  distribution  studies  were  carried  out  with  instrumentation  mounted 
in  the  bottom  window.  The  top  window  was  used  as  an  observation  port  for  photography  of  oil  flow  patterns 

2.3  Surface  Measurements 


Oil  flow  studies:  A  technique  developed  for  the  present  test  conditions  used  a  fluorescent  oil 
applied  upstream  of  the  region  of  study  before  the  tunnel  was  started.  Approximately  10  seconds  after  the 
tunnel  was  operating,  photographs  were  taken  of  the  oil  pattern.  The  particular  mixture  used  consisted  of 
clear  motor  oil  and  green  fluorescent  dye  thinned  with  kerosene.  The  oil  patterns  were  illuminated  by  a 
high  intensity  ultra  violet  light  source.  This  particular  combination  of  ultra  violet  excitation  and 
green  fluorescence  made  it  possible  to  filter  out  the  ultra  violet  reflections  coming  from  the  metal  test 
surface.  Good  results  were  obtained  under  all  test  conditions. 

Surface  static  pressure:  Measurements  were  made  on  Model  1  using  an  instrumented  plate  which  fitted 
into  the  bottom  window  port.  A  large  number  of  static  pressure  taps  was  necessary  to  map  the  full  flow 
field  involving  the  interaction  and  to  account  for  the  change  in  shock  location  when  the  shock  generator 
angle  was  varied.  During  a  test,  a  fixed  shock  generator  position  was  held  while  a  scanivalve  measured 
the  static  pressure  distribution.  The  basic  static  pressure  window  port  could  not  be  used  for  Model  2,  so 
the  test  surface  itself  was  instrumented  with  an  arrangement  of  static  pressure  orifices. 

Heat  transfer:  Several  techniques  were  investigated  to  measure  the  heat  transfer  rate  over  the  test 
area.  The  method  finally  chosen  was  a  quasi-transient  method  using  a  slug  calorimeter  since  it  matched 
best  to  the  tunnel  operating  procedures,  gave  reasonable  spatial  resolution  and  is  rather  simple  to  manu¬ 
facture.  Figure  3  shows  a  photograph  of  the  heat  transfer  plug  and  the  insert  which  provides  the  hot  air. 
A  line  drawing  of  the  slug  calorimeter  arrangement  is  given  in  Fig.  4. 

To  determine  the  heat  transfer  rate  from  the  wall  to  the  gas  at  uniform  wall  temperature  conditions 
it  is  necessary  that  the  test  piece  (slug)  be  heated  to  a  temperature  higher  than  the  surrounding  wall 
temperature.  This  heating  was  provided  by  impinging  a  tiny  jet  of  hot  air  to  the  inside  of  each  slug 
calorimeter.  After  the  hot  air  heating  is  stopped,  the  slug  would  attempt  to  reach  the  wall  recovery 
temperature  as  time  progresses.  At  an  intermediate  time,  t0  ,  while  the  slug  temperature  is  dropping  it 
equals  to  the  wall  temperature.  From  an  exponential  curve  fitted  to  the  measured  slug  temperature-time 
data  points,  the  temperature- time  gradient  was  evaluated  at  t0  .  The  convective  heat  transfer  follows 
from  this  temperature-time  derivative  and  the  calorimeter  properties.  Heat  conduction  from  the  slug  to 
the  surrounding  wall  was  minimized  by  separating  the  slugs  from  the  aluminum  plug  by  nylon  insulation.  It 
was  determined  that  no  correction  due  to  insulation  losses  were  necessary. 

Skin  friction:  The  Preston  tube  was  used  to  get  some  indication  of  the  skin  friction.  As  the  pre¬ 
sent  problem  involves  both  pressure  gradients  and  cross  flows,  there  will  be  uncertainties  in  the  results, 
particularly  for  the  case  of  large  cross  flows.  The  Preston  probe  was  first  aligned  with  the  surface  oil 
flow  direction,  then  rotated  10°  in  both  directions.  Although  the  Preston  tube  reading  was  relatively 
insensitive  to  such  rotation,  the  largest  measured  value  was  used  to  evaluate  the  skin  friction.  The  cali¬ 
bration  scheme  proposed  by  Bradshaw7  was  used  to  evaluate  the  skin  friction.  He  extended  the  subsonic 
calibration  equation  to  compressible  flow  in  an  alternate  way.  Instead  of  using  calibration  factors  to 
transform  to  the  subsonic  calibration  coordinates,  he  added  a  compressibility  correction  function  to  the 
subsonic  calibration  equation.  This  new  function  has  variables  which  are  based  upon  wall  conditions  only. 

2 . 4  Flow  Field  Measurements 


To  carry  out  detailed  flow  field  measurements  in  the  interaction  region,  a  probe  drive  system  with 
two  degrees  of  freedom  was  used.  It  moved  the  probes  normal  to  the  wall  and  yawed  them.  The  probes  them¬ 
selves  were  supported  by  a  streamlined  section  which  extended  from  the  top  wall  down  towards  the  bottom 
wall  boundary  layer  where  the  measurements  were  taken. 

Pitot-yaw  surveys:  To  determine  the  flow  angles  in  the  X-Y  plane  throughout  the  flow  field,  a  cobra 
Pr°be  has  been  employed.  The  cobra  probe  was  traversed  and  yawed  in  such  a  way  that  the  two  side  pressures 
remained  equal  so  that  the  flow  yaw  angle  and  pitot  pressure  could  be  determined  at  the  same  time. 


Static  pressure  surveys:  Since  it  was  clear  that  significant  variations  in  the  static  pressure  would 
take  place  within  the  interaction  region,  it  was  necessary  to  determine  the  local  static  pressure  in  order 
to  evaluate  the  local  flow  velocities  and  density.  To  obtain  a  measure  of  the  local  values,  a  small  cone- 
cylinder  probe  was  designed  to  be  used  in  conjunction  with  the  probe  drive  mechanism  noted  previously. 

The  static  pressure  orifices  were  located  11  diameters  downstream  of  the  tip.  The  probe  rotated  in  the 
X-Y  plane  about  the  orifice  location  and  was  kept  aligned  with  the  known  flow  direction.  A  calibration 
of  the  angle  sensitivity  of  the  probe  in  uniform  flow  was  performed  and  it  was  found  that  over  *  4°,  the 
static  pressure  varied  less  than  2%  from  the  correct  value. 

Total  temperature  surveys:  To  complete  the  measurements  from  which  the  local  velocities  can  be  ob¬ 
tained,  measurements  were  made  of  the  total  temperature  throughout  the  interaction  region.  Although  the 
wall  temperature  was  approximately  15%  above  its  adiabatic  value,  the  measurement  of  small  variations  of 
total  temperature  provide  for  increased  accuracy  of  the  deduced  flow  field  data.  The  probe  used  was  a 
fine  wire  probe  which  had  a  thermocouple  junction  in  the  middle  of  the  wire.  This  thermocouple  consisted 
of  the  junction  between  alumel  and  chromel  wires  which  had  a  thickness  of  0.002  inches  diameter. 

A  more  detailed  description  of  the  experimental  techniques  used  is  given  in  Ref.  5. 

3.  INCOMING  BOUNDARY  LAYER  CHARACTERISTICS 

The  boundary  layers  that  interacted  with  the  shock  wave  were  surveyed  extensively  and  were  shown  to 
be  two  dimensional  across  the  width  of  the  wind  tunnel.  It  should  be  noted,  however,  that  the  upstream 
boundary  of  the  interaction  region,  as  determined  from  the  start  of  the  wall  static  pressure  increase,  is 
a  swept  line  which  is  approximately  parallel  to  the  shock  wave.  This  results  in  a  variation  in  thickness 
of  the  boundary  layer  coming  into  the  interaction  region  since  the  boundary  layer  is  growing  in  the  stream- 
wise  direction.  The  boundary  layer  that  interacted  with  the  shock  wave  in  the  Model  1  configuration  ori¬ 
ginated  from  the  tunnel  nozzle  and  increased  in  thickness  from  0.52  to  0.58  inches  along  a  typical  inter¬ 
action  region  for  a  10  degree  shock  generator  angle.  The  boundary  layer  for  Model  2  was  generated  on  a 
flat  test  plate  which  spanned  the  tunnel  and  was  located  2  inches  above  the  tunnel  bottom  wall.  The  cor¬ 
responding  increase  in  thickness  of  this  thinner  Model  2  boundary  layer  was  much  larger  because  of  its 
shorter  running  length  and  amounted  to  a  variation  of  0.11  to  0.16  inches.  These  increases  of  the  incoming 
boundary  layer  thickness  along  the  shock  wave  have  to  be  taken  into  account  when  the  scale  of  the  inter¬ 
action  length  is  evaluated  in  terms  of  a  local  boundary  layer  thickness. 

Boundary  layer  parameters  derived  from  the  measurements  are  tabulated  for  two  typical  initial  pro¬ 
files  for  Model  1  and  2,  respectively: 


6  (inches) 

6* (inches) 

8 (inches) 

cf 

Reinch 

Model  1 

0.540 

0.162 

0.0S2 

0.00119 

1.6  x  106 

Model  2 

0.132 

0.040 

0.008 

0.00150 

1.6  x  106 

The  corresponding  velocity  profiles  were  transformed  to  the  incompressible  law  of  the  wall  variables  by 
using  the  van  Driest  transformation  as  given  in  Ref.  8.  These  profiles  are  then  compared  with  the  incom¬ 
pressible  wall-wake  law  as  given  by 


V*Ftn(f)  +  C-iw(£) 


(1) 


where 


N(£)  «  1  -  cos(r  £),  C  =  5.0,  k  =  0.41,  n  =  0.52  . 


The  comparison  (see  Fig.  5)  shows  that  the  test  boundary  layers  can  be  described  as  fully  developed, 
"equilibrium",  turbulent  boundary  layers. 


The  overall  size  of  the  region  of  study  for  both  Model  1  and  2  in  terms  of  boundary  layer  thickness 
is  represented  in  Fig.  6  and  compared  with  those  of  previous  studies  by  McCabe,  Lowrie  and  Token.  Model  1 
provides  results  in  an  area  which  is  also  covered  by  McCabe  and  Lowrie,  but  outside  the  region  studied  by 
Token.  Model  2  covers  a  larger  range  not  heretofore  explored. 


The  general  coordinate  system  used  in  obtaining  and  presenting  the  data  is  also  shown  in  Fig.  6. 


The  shock  location,  calculated  from  the  oblique  shock  theory,  is  used  as  the  origin  for  the  X  coordinate. 


4.  ANALYSIS  OF  RESULTS  AND  DISCUSSION 


4. 1  Three-Dimensional  Separation 


The  notion  of  flow  separation  on  planar  surfaces  in  both  two-dimensional  as  well  as  three-dimensional 
flow  is  that  a  boundary  layer,  flowing  tangentially  alongside  a  solid  surface,  breaks  away  from  this  sur¬ 
face  due  to  an  adverse  pressure  gradient.  In  planar,  two-dimensional  flows  the  customary  definition  of 
separation  requires  that  the  mean  wall  shear  stress  vanishes  at  the  separation  point.  This  definition  is 
quite  satisfactory  for  such  cases  as  it  implies  the  start  of  flow  reversal.  In  three-dimensional  flow, 
however,  the  situation  is  more  complicated,  illustrated  by  the  fact  that  three-dimensional  separation  has 
a  two-fold  mechanism  (see  Lighthill,  Ref.  9).  This  mechanism  is  explained  by  using  the  concept  of  limiting 
surface  streamlines  (or  skin  friction  lines).  If  one  considers  a  streamline  at  a  small  distance  from  the 
surface,  then  it  follows  from  the  continuity  equation  that  this  streamline  has  to  depart  from  the  surface 
not  only  in  case  of  vanishing  wall  shear-stress,  but  also  in  a  case  where  the  topography  of  the  limiting 


surface  streamlines  is  such  that  they  converge.  They  can  coalesce  and  form  a  line  of  ordinary  separation. 
Although  the  notion  of  flow  separation  as  applied  to  two-  or  three-dimensional  flow  is  the  same,  the  main 
mechanism  through  which  flow  separation  takes  place  can  be  different;  for  example,  three-dimensional  flow 
can  separate  from  the  wall  along  a  separation  line  without  the  mean  surface  shear  stress  becoming  zero. 

This  type  of  flow  separation  is  known  as  ordinary  separation  (see  Maskell,  Ref.  10)  and  has  to  be  dis¬ 
tinguished  from  singular  separation  whose  main  characteristic  is  that  the  mean  wall  shear  stress  vanishes 
at  the  separation  point. 

In  applying  these  ideas  of  flow  separation  to  the  present  problem,  previous  investigators  have  pro¬ 
posed  several  ways  in  which  the  onset  of  flow  separation,  over  any  region  of  significant  size,  can  be 
determined  from  surface  flow  patterns.  The  main  characteristic  of  a  line  of  ordinary  separation  is  the 
tangential  convergence  of  the  surface  streamlines  into  one  single  line.  Maskell  proposed  that  the  order 
of  magnitude  of  the  rate  of  this  convergence  be  the  same  on  each  side  of  the  separation  line.  Rogers  and 
Hall^*  have  relaxed  this  criterian  for  separation,  by  suggesting  that  for  the  condition  of  incipient  sepa¬ 
ration  the  limiting  surface  streamlines  converge  toward  the  separation  line  only  in  the  upstream  region, 
immediately  downstream  they  are  parallel  to  it  and  further  downstream  they  turn  progressively  away  from  it. 

Although  these  criteria  of  ordinary  flow  separation  are  based  upon  limiting  surface  streamlines,  and 
thus  have  to  be  verified  by  actual  flow  field  measurements,  they  are  believed  to  be  physically  correct. 
Care,  however,  must  be  taken  in  applying  these  criteria  especially  in  turbulent  boundary  layers  where  the 
equivalence  between  limiting  surface  streamlines  (skin  friction  lines)  and  oil  flow  lines  is  not  an  obvious 
matter.  In  the  following  sections,  there  are  several  discussions  specifically  aimed  at  trying  to  clarify 
questions  about  this  ordinary  flow  separation.  Finally,  objections  will  be  raised  against  McCabe's  cri¬ 
terion  for  separation  as  used  in  other  investigations  (Refs.  1-4)  of  the  present  problem. 

4.2  Surface  Flow  Patterns 

Some  results  obtained  during  the  oil  flow  studies  are  presented  in  Figs.  7  through  10  with  ,G  =  4 
and  10  degrees  for  both  Model  1  and  2.  The  oil  lines  for  aG  =  4°  downstream  of  the  calculated  shock 
location  are  deflected  to  a  maximum  angle  of  about  10°,  which  is  about  twice  the  shock  generator  angle 
but  remains  well  below  the  shock  wave  angle.  For  aG  =  10°  the  oil  lines  turn  to  angles  which  are  sub¬ 
stantially  larger  than  the  shock  wave  angle  and  coalescence  of  the  oil  streaks  occurs  as  a  consequence. 
Special  attention  should  be  paid  to  the  difference  in  scale  of  an  average  boundary  layer  thickness  between 
Model  1  and  2.  Because  of  this  large  region  of  Model  2  (see  Fig.  6)  there  emerges  an  area  downstream  of 
the  shock  wave  for  aG  =  10°  where  the  oil  streak  lines  are  parallel  to  the  shock  generator  again.  This 
region  (II)  is  indicated  on  Fig.  10. 

The  interpretation  is  that  the  interaction  is  completed  and  the  boundary  layer  has  again  reached  a 
planar,  two-dimensional  character  in  this  area  indicated.  The  implication  of  this  observation  is  that, 
for  example,  in  the  case  of  aG  *  10°,  Y  «  2.2S  inches,  it  takes  about  30  boundary  layer  thicknesses  down¬ 
stream  of  the  shock  wave  before  all  cross  flows  have  disappeared  from  the  boundary  layer.  This  length  of 
the  downstream  extent  of  the  interaction  has  never  been  observed  in  previous  investigations  and  is  a  direct 
result  of  the  large  region  covered  by  Model  2.  The  fact  that  the  end  of  the  interaction  does  not  appear  on 
the  surface  flow  patterns  for  oG  smaller  than  6°  is  an  indication  that  the  downstream  extent  of  the 
interaction  in  the  X-direction  is  larger  for  smaller  'hock  generator  angles. 

These  surface  flow  patterns  can  be  reduced  to  quantitative  data  by  measuring  the  angles  of  local  oil 
lines  with  respect  to  the  X-direction  along  a  line  of  Y  =  constant.  Results  obtained  from  that  process 
for  shock  generator  angles  ranging  from  2  up  to  12  degrees  are  presented  in  Fig.  11  for  Model  2.  The 
theoretical  shock  wave  angles  are  shown  to  the  left  hand  side  of  these  figures.  From  Fig.  11  it  can  be 
seen  that  the  rate  of  decrease  of  the  local  oil  line  angles  with  distance  Xs  for  small  generator  angles, 
xG  <_  6°,  is  much  smaller  than  for  the  larger  aG's  •  This  results  in  a  larger  downstream  extent  for 
weaker  shocks.  The  explanation  of  this  seemingly  contradictory  result  lies  in  the  topography  of  the  flow 
field  that  is  associated  with  the  surface  flow  patterns. 

For  xG  =  6°  and  smaller,  the  total  yaw  angles  occurring  in  the  flow  field  remain  smaller  than  the 
shock  wave  angle.  These  flow  fields  are  similar  to  a  classical  three  dimensional  boundary  layer  flow  that 
negotiates  a  transverse  pressure  gradient.  This  transverse  pressure  gradient  causes  the  slower  moving 
fluid  in  the  bottom  portion  of  the  boundary  layer  to  deflect  to  larger  angles  than  the  faster  moving  fluid 
in  the  outer  portion  of  the  boundary  layer.  This  process  in  which  cross  flows  are  produced  by  the  trans¬ 
verse  pressure  gradient  is  an  instantaneous  process  or  in  other  words  the  cross  flows  are  created  at  the 
physical  location  where  the  transverse  pressure  gradients  are  present.  These  pressure  gradients  disappear 
in  going  downstream  of  the  shock  location.  The  cross  flows,  however,  do  not  vanish  instantaneously,  but 
rather  decay  slowly.  The  driving  force  in  this  decay  process  is  the  magnitude  of  the  cross  flow  itself. 

So  as  the  cross  flow  gets  smaller,  the  rate  of  decay  decreases,  and  the  downstream  limit  of  no  cross  flow 
is  only  reached  asymptotically,  resulting  in  very  large  downstream  extent  for  aG  <_  6°. 

For  larger  generator  angles  (aG  8°)  other  processes  play  a  role  in  the  decay  of  cross  flows. 

Since  the  yaw  angles  are  no  longer  smaller  than  the  shock  angle,  it  can  no  longer  be  considered  as  regular 
flow  field,  but  rather  a  flow  with  large  secondary  flows  imbedded  in  it.  The  fluid  in  which  the  large 
cross  flows  are  produced  is  transported  along  the  shock  direction  and  thus  out  of  the  plane  of  observation, 
the  X-Z  plane.  This  reorganization  process  of  the  flow,  as  described  in  a  later  section  on  the  flow  model, 

causes  fluid  in  which  much  smaller  cross  flows  are  present  to  pitch  down  toward  the  wall  resulting  in  a 

strong  decay  of  oil  line  angles  with  distance  Xs  for  aG  ^8°. 

McCabe1  has  proposed  a  simple  approximate  secondary  flow  theory.  This  theory  assumes  that  the 
limiting  streamline  angle  does  not  vary  with  distance  downstream  of  the  shock  wave.  This  is  an  oversimpli¬ 
fication  of  the  actual  situation  as  can  be  seen  from  Fig.  11.  However,  if  the  predicted  limiting  stream¬ 
line  angle  is  assumed  to  be  a  representation  of  the  maximum  of  the  distribution  of  oil  line  angles  for  a 
given  shock  generator  deflection,  aojiax  ’  then  the  theory  is  a  good  approximation  of  the  present  experi¬ 
mental  data,  as  shown  in  Fig.  12.  The  theoretical  shock  wave  angle  is  also  indicated  in  the  figure. 


McCabe  has  also  suggested  a  criterion  for  ordinary  separation  which  defines  the  flow  as  separated  if 
the  oil  lines  are  deflected  to  an  angle  larger  than  the  shock  wave  angle.  This  criterion  indicates  that 
the  flow  is  separated  for  a q  >  7.S0  from  the  present  data  (see  Fig.  12).  The  justification  of  this  cri¬ 
terion,  however,  is  very  much  in  question  since  it  would  not  coincide  with  the  first  appearance  of  an 
ordinary  separation  line  as  shown  by  McCabe  himself.  This  criterion  also  implies  that  the  minimum  shock 
generator  angle  needed  for  flow  separation  decreases  indefinitely  as  the  shock  wave  angle  decreases  with 
increasing  Mach  number  (see  Korkegi,  Ref.  12).  This  last  observation  indicates  that  even  the  practical 
value  of  McCabe's  critericn  is  not  at  all  clear. 

4. 3  Static  Pressure  Distributions 

A  typical  static  pressure  distribution  is  shown  in  Fig.  13.  The  distance  Xs  is  normalized  by  the 
boundary  layer  thickness,  Sj  ,  at  the  beginning  of  the  pressure  rise.  The  upstream  extent  of  the  inter¬ 
action  region  measured  from  the  shock  position  is  about  12  6j  and  largely  independent  of  the  shock  strength 
at  this  station  of  Y  «  2.25  inches  of  Model  2.  The  downstream  extent  is  more  difficult  to  define  in  terms 
of  the  static  pressure  distribution  because  this  limit  is  reached  asymptotically.  The  observation  that 
the  scale  of  this  downstream  extent  is  of  the  order  of  20  to  30  times  the  incoming  boundary  layer  thickness 
at  Y  =  2.25  inches  is  of  quite  some  significance.  This  scale  was  also  extracted  from  the  surface  flow 
patterns . 

Static  pressure  distributions  obtained  along  several  instrumentation  lines  of  Y  =  constant  were  also 
cross  plotted  as  isobar  patterns.  A  typical  isobar  pattern  is  shown  in  Fig.  14  for  ag  =  10°,  Model  1. 

As  can  be  noted,  the  lines  of  constant  static  pressure  are  approximately  parallel  to  the  shock  wave  and 
only  to  a  small  amount  diverging  along  the  shock  direction.  This  indicates  that  the  pressure  distribution 
is  quasi -two- dimensional  in  the  sense  that  this  flow  quantity  is  constant  in  the  shock  direction  to  a  good 
approximation. 

Some  comparison  can  be  made  with  the  data  obtained  by  other  investigators.  The  test  conditions  are 
listed  in  the  following  table: 


Mach 

Number 

6 

ave 

(inches) 

Re.  . 
inch 

Re6 

ave 

Lowrie 

(Ref. 2) 

3.44 

0.412 

1.02  x  106 

4.20  x  105 

McCabe 

(Ref.  1) 

2.94 

0.230 

2.06  x  105 

0.47  x  105 

Token 

(Ref. 3) 

3.71 

6.0 

2.91  x  105 

17.5  x  105 

Model 

1 

2.95 

0.55 

1.6  x  106 

8.8  x  105 

Model 

2 

2.95 

0.133 

1.6  X  106 

2.1  x  105 

The  test  Mach  number  of  all  these  studies  is  around  3  and  the  wall  temperatures  are  all  near  adia¬ 
batic.  The  boundary  layer  thickness,  however,  varies  by  a  factor  of  40.  If  the  intersection  of  a  line 
drawn  tangent  to  the  maximum  slope  of  a  static  pressure  distribution  and  the  line  P  *  1  ,  is  taken  as  the 
beginning  of  interaction,  then  the  upstream  extent,  AX  ,  can  be  defined  as  the  distance  between  the  calcu¬ 
lated  shock  location  and  this  beginning.  The  upstream  extent  of  the  interaction,  nondimensionalized  by 
the  local  incoming  boundary  layer  thickness,  is  shown  in  Fig.  15.  The  Y  coordinate  in  this  plot  is  non¬ 
dimensionalized  by  the  boundary  layer  thickness,  6j  ,  averaged  over  the  distance  between  the  shock  gener¬ 
ator  leading  edge  and  the  station  Y  -  constant  considered.  This  average  of  6j  is  indicated  by  6avg  . 

The  general  agreement  of  the  data  in  Fig.  15,  despite  the  difference  in  Mach  number  and  deflection 
angles,  has  an  important  implication.  If  one  recalls  that  6j  varies  by  a  factor  of  40,  one  can  reach 
the  conclusion  that  the  boundary  layer  thickness  can  provide  a  basis  for  scaling  parameters  in  this  type 
of  problem. 

4.4  Heat  Transfer  Distribution 

The  heat  transfer  results  are  presented  as  a  nondimensional  ratio  of  the  local  heat  transfer  coeffi¬ 
cient  divided  by  the  flat  plate  value  predicted  by  the  van  Driest  method^  for  the  test  conditions  upstream 
of  the  interaction.  It  should  be  noted  that  the  local  heat  transfer  coefficient  is  defined  as  the  local 
heat  transfer  rate  nondimensionalized  by  free  stream  conditions  ahead  of  the  interaction.  The  wall  temper¬ 
ature  conditions  are  uniform  throughout  the  interaction  region. 

The  results  indicate,  as  shown  in  Fig.  16,  that  heat  transfer  rate  decreases  somewhat  as  the  shock 
is  approached  in  the  streamwise  direction,  and  increases  almost  linearly  with  the  downstream  distance  from 
the  shock  position.  No  conclusion  about  the  peak  heat  transfer  could  be  reached  as  the  maximum  occurred 
at  the  boundary  of  the  region  of  study  of  Model  1.  The  heat  transfer  results  are  also  cross  plotted  as  a 
heat  transfer  pattern  in  the  X-Y  plane  and  an  example  is  given  in  Fig.  17  for  aG  »  10°,  Model  1.  The 
results  show  a  general  uniformity  of  the  heat  transfer  field  and  demonstrate  again  the  quasi-two-dimensional 
nature  of  the  interaction  region. 

4 . 5  Preston  Tube  Measurements 

The  skin  friction  was  deduced  from  Preston  tube  measurements  taken  along  Y  =  4.0  inches  on  Model  1 
for  aG  «  2°  to  10°.  These  results  of  these  measurements  are  shown  in  Fig.  18  where  the  skin  friction  has 
been  nondimensionalized  by  the  upstream  flat  plate  value.  The  most  significant  observation  is  that  the 
skin  friction  distribution  is  almost  equivalent  to  the  heat  transfer  distribution.  This  implies  that  the 
Reynolds  analogy  is  valid  as  a  first  approximation.  If  ordinary  separation  occurs  it  is  not  required  that 
the  skin  friction  vanish,  therefore,  no  conclusion  about  separation  can  be  reached  by  these  results. 


40- 


4.6  Analysis  of  Complete  Flow  Field  Data 

A  complete  picture  of  the  mean  flow  field  in  the  X-Z  plane  can  be  constructed  from  the  measured 
quantities,  pitot  pressure,  vaw  angle,  static  pressure,  and  total  temperature.  If  the  static  pressure 
and  the  pitot  pressure  at  a  point  in  the  field  are  known,  the  Mach  number  can  be  calculated  directly  from 
the  isentropic  relations  for  subsonic  conditions,  and  from  the  Rayleigh  Pitot  formula  in  an  iterative 
cycle  for  supersonic  cases.  The  total  energy  equation  was  then  applied  to  give  the  static  temperature  and 
thus  the  density  and  u,v  velocity  components.  In  addition,  the  velocity  component  in  the  Z-direction, 
w  ,  can  be  deduced  from  the  density  and  u,v  velocity  distributions  by  integrating  the  continuity  equa¬ 
tion  and  using  the  boundary  condition  w  =  0  at  4=0.  Hence' 


w  =  -  £  fo  {h  (DU)  +  w  (pv)!  dz 


(2) 


The  derivatives  in  the  integrand  were  obtained  by  locally  fitting  a  second  order  polynomial  through  three 
data  poi..ts  followed  by  a  differentiation  of  the  polynomial.  Data  was  taken  in  one  X-Z  plane  only.  To 
overcome  this  difficulty  of  evaluating  the  derivative,  the  data  in  the  X-Z  plane  at  Y  =  4.0  inches 

3y 

was  projected  along  the  shock  direction  to  higher  and  lower  values  of  Y  .  This  procedure  assumes  that 
the  interaction  region  is  perfectly  two-dimensional  along  the  shock  direction,  but  it  was  found  that  the 
resultant  w  was  not  sensitive  to  small  changes  of  the  projection  direction,  e.g.,  a  slight  divergence 
of  the  projection  directions  through  the  interaction.  The  accuracy  of  the  w  component  obtained  by  dif¬ 
ferentiation  and  integration  of  the  original  data  is  less  than  that  of  u  and  v  .  After  the  w-component 

w 

is  computed, the  pitch  angle  8  (defined  as  8  =  arctan  (— ))  can  be  found.  It  should  be  noted  that  the 

u  2  2 

measurement  technique  of  obtaining  u  and  v  is  only  valid  if  w  remains  small  compared  to  ’’u'-  +  v 
This  requirement  was  checked  a  posteriori  and  found  to  be  satisfied. 

The  static  pressure  distribution  in  the  X-Z  plane  at  Y  =  4.0  inches  for  ng  =  10°,  shown  in  Fig.  19, 
indicates  that  within  the  interaction  there  is  no  major  region  (in  terms  of  6)  where  the  static  pressure 
is  uniform.  The  flow  ahead  of  the  shock  is  characterized  by  compression  waves  extending  to  well  outside 
the  incoming  boundary  layer.  The  static  pressure  field  downstream  of  the  shock  has  a  more  complicated 
distribution.  The  Mach  number  distribution  in  this  X-Z  plane  is  mostly  supersonic,  Fig.  20.  A  shallow 
subsonic  region  was  measured  close  to  the  wall  near  the  shock  location. 

Large  yaw  angles  are  measured  in  the  interaction  region  (see  Fig.  21).  The  largest  yaw  angles  of 
55°  occurred  a  small  distance  from  the  wall  and  Xs  <  3  .  Despite  this  small  region  where  the  yaw  angles 
are  larger  than  the  shock  wave  angle,  there  is  an  extensive  region  downstream  of  the  shock  wave  where  the 
yaw  angles  are  substantially  less  than  the  shock  generator  angle  (ag  =  10°). 

The  pitch  angle  8  through  the  interaction  (shown  in  Fig.  22)  provides  some  additional  information. 
The  compression  waves  ahead  of  the  shock  do  not  only  cause  increasing  yaw  angles  but  also  increasing  pitch 
angles.  Going  through  the  shock  wave  these  pitch  angles  don’t  seem  to  decrease  initially,  however,  at 
Xs  3  rapid  changes  take  place  and  the  pitch  angle  becomes  negative. 

In  view  of  the  question  about  ordinary  separation,  it  should  be  noted  from  this  figure  that  the  pitch 
angle  tends  to  zero  as  the  wall  is  approached.  This  observation  is  of  particular  importance  at  the  posi¬ 
tion  (Xs  \  -3)  where  tne  oil  lines  tend  to  coalesce.  This  is  the  region  where  one  would  expect  locally 
large  pitch  angles.  The  surveys  were  taken  at  stations  one  boundary  layer  thickness,  6j  ,  apart.  Keeping 
this  streamwise  resolution  in  mind  it  can  be  concluded  that  a  local  coalescence  of  oil  flow  lines  is  not 
necessarily  related  to  locally  large  streamline  inclination  with  respect  to  the  wall  surface.  There  is, 
however,  a  tendency  of  the  flow  to  depart  from  the  wall  (8  positive)  but  it  is  uniformly  spread  out  over 
many  boundary  layer  thicknesses. 

4.7  Flow  Models  and  Separation 

To  depict  the  flow  field  in  the  X-Z  plane  at  Y  =  4.0  inches,  the  following  flow  model  is  constructed 
for  xg  =  10°,  see  Fig.  23.  The  region  A  is  the  undisturbed  uniform  flow  ahead  of  the  interaction  with  a 
planar  turbulent  boundary  layer  at  the  wall.  The  typical  character  of  region  B  is  that  it  has  very  slowly 
increasing  yaw  angles  which  remain  lower  than  the  shock  generator  angle  (ag  =  10°)  even  in  subregions  Bx 
and  Bj  and  only  reach  the  10°  level  far  downstream.  In  region  C  (the  area  extending  downward  from  region 
B  and  approaching  the  wall)  large  yaw  angles  prevail  (10°  ^  a  55°) .  Yaw  angles  larger  than  the  shock 
wave  angle  (xsf,  =  28°)  occur  in  subregion  Cj.  The  far  flow  field  downstream  of  the  shock  and  outside 
region  B  contains  waves  originating  from  the  interaction,  region  D.  The  decrease  in  pitch  angles  2  ob¬ 
served  by  going  from  subregion  Bj  to  Bj  generates  expansion  waves  in  the  external  flow  (subregion  D j ) . 
Compression  waves  are  present  in  subregion  Dj  where  the  flow  has  to  turn  parallel  to  the  wall  again. 

The  corresponding  flow  model  for  ag  =  4°  is  also  presented,  see  Fig.  24.  The  measurements  on  which 
this  model  is  based  are  given  in  Ref.  5.  Unlike  the  ag  =  10°  case,  the  yaw  angles  in  this  case  remain 
well  below  the  shock  wave  angle.  This  means  that  no  reorganization  process  occurs  as  in  the  case  of  1(-  = 
10°  where  the  net  transport  of  fluid  along  the  shock  direction  incurred  negative  values  of  8  .  Although 
it  was  determined  that  the  maximum  pitch  angles  were  quite  small  (1.5°),  they  are  significantly  larger  than 
those  associated  with  boundary  layers  developing  under  zero  pressure  gradient  conditions. 

It  is  obvious  that  for  ag  =  4°,  although  the  streamlines  are  slightly  departing  from  the  wall 
(S  <_  1.5°),  that  there  is  no  question  about  separation.  For  the  xg  =  10°  case,  however,  the  question  of 
whether  the  boundary  layer  should  be  called  separated  is  not  an  obvious  matter.  If  one  examines  the  oil 
flow  patterns  for  ag  *  10°  Figs.  8  and  10)  more  carefully,  the  observation  can  be  made  that  although  the 
oil  lines  are  converging  and  tend  to  coalesce  upstream  of  the  calculated  shock  position,  they  do  not  con¬ 
verge  into  one  single  line.  The  region  in  which  coalescence  occurs  is  growing  in  size  in  the  direction 
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along  the  shock.  This  observation  is  consistent  with  the  idea  that  the  site  of  the  interaction  region  is 
increasing  with  distance  along  the  shock.  It  is  exactly  this  gradual  increase  in  site  that  permits  the 
oil  lines  to  converge  without  incurring  locally  large  pitch  angles  B  corresponding  to  separation  of  the 
flow  from  the  wall.  This  directly  leads  to  the  conclusion  that  the  flow  is  also  unseparated  for  ir  =  10°. 

The  results  clearly  indicate  that  the  qualitative  difference  between  attached  flows  and  flows  with 
ordinary  separation  for  the  present  problem  cannot  be  taken  for  granted.  The  definition  of  ordinary  sepa¬ 
ration  becomes  then  a  matter  of  some  conjecture  but  does  not  seem  to  be  as  crucial  as  its  counterpart  in 
planar  two-dimensional  flow  concerning  singular  separation. 

5.  CONCLUSIONS 

A  detailed  study  was  carried  out  of  the  three-dimensional  flow  fields  caused  by  the  interaction  of 
oblique  shock  waves  and  a  planar  turbulent  boundary  layer.  The  main  results  are  the  following: 

1.  Two  flow  field  models  have  been  constructed  based  on  experimental  data  for  shock  generator 
angles  of  4°  and  10°.  The  measurements  give  a  complete  description  of  the  flow  fields. 

2.  Based  on  the  10°  flow  field  data,  it  is  concluded  that  McCabe's  criterion  for  incipient  separa¬ 
tion  is  not  sufficient  to  define  flow  separation. 

3.  The  interaction  region  is  quasi-two-dimensional  in  an  area  at  sufficient  distance  from  the  shock 
generator.  This  means  that  although  all  three  velocity  components  play  an  important  role,  the  fluid  vari¬ 
ables  vary  only  slowly  along  the  shock  direction. 

4.  The  heat  transfer  rate  decreases  as  the  shock  is  approached  in  the  streamwise  direction  and 
increases  almost  linearly  with  downstream  distance  from  the  shock  location. 

5.  Reynolds  analogy  is  valid  as  a  first  approximation  to  the  problem. 

6.  Boundary  layer  thicknesses  will  provide  a  basis  for  scaling  parameters  of  the  upstream  extent 
as  a  function  of  the  transverse  distance  Y. 

7.  The  downstream  extent  of  the  interaction  in  X-direction  is  at  least  four  times  as  large  as  the 
upstream  one. 
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Fig.  16:  Heat  Transfer,  Y  =  4.0  inches.  Model  1 


Fig.  17:  Heat  Transfer  Pattern, 
c>G  =  10°,  Model  1 


Fig.  18:  Skin  Friction, 

Y  =  4.0  inches,  Model  1 
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Fig.  19:  Static  Pressure 
Contour  Plot  in  XZ- 
Plane,  Y  =  4.0  inches, 
aG  =  10°,  Model  1 
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Fig.  20:  Mach  Number 
Contour  Plot  in  XZ- 
Plane,  Y  =  4.0  inches 
^  =  10°.  Model  1 


